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Abstract 

Wightman function, the vacuum expectation values of the field square and the energy- 
momentum tensor are investigated for a massive scalar field with an arbitrary curvature 
coupling parameter in the region between two infinite parallel plates moving by uniform 
proper acceleration. We assume that the field is prepared in the Fulling-Rindler vacuum 
state and satisfies Robin boundary conditions on the plates. The mode-summation method 
is used with a combination of a variant of the generalized Abel-Plana formula. This allows 
to extract manifestly the contributions to the expectation values due to a single boundary 
and to present the second plate- induced parts in terms of exponentially convergent integrals. 
Various limiting cases are investigated. The vacuum forces acting on the boundaries are 
presented as a sum of the self-action and 'interaction' terms. The first one contains well 
known surface divergences and needs a further renormalization. The 'interaction' forces 
between the plates are investigated as functions of the proper accelerations and coefficients 
in the boundary conditions. We show that there is a region in the space of these parameters 
in which the 'interaction' forces are repulsive for small distances and attractive for large 
distances. 



PACS number(s): 03.70.+k, ll.10.Kk, 04.62,+v 



1 Introduction 

The Casimir effect is one of the most interesting macroscopic manifestations of the nontrivial 
structure of the vacuum state in quantum field theory (see, e.g., P3C3IS1IIJ and references therein). 
The effect is a phenomenon common to all systems characterized by fluctuating quantities and 
results from changes in the vacuum fluctuations of a quantum field that occur because of the 
imposition of boundary conditions or the choice of topology. It may have important implications 
on all scales, from cosmological to subnuclear, and has become in recent decades an increasingly 
popular topic in quantum field theory. It is well known that the uniqueness of the vacuum state is 
lost when we work within the framework of quantum field theory in a general curved spacetime or 
in non-inertial frames. In particular, the use of general coordinate transformations in quantum 
field theory in flat spacetime leads to an infinite number of unitary inequivalent representations 
of the commutation relations. Different inequivalent representations will in general give rise to 
different vacuum states. For instance, the vacuum state for a uniformly accelerated observer, the 
Fulling-Rindler vacuum 0013111121) turns out to be inequivalent to that for an inertial observer, 
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the familiar Minkowski vacuum. Quantum field theory in accelerated systems contains many 
special features produced by a gravitational field. In particular, the near horizon geometry 
of most black holes is well approximated by Rindler spacetime and a better understanding 
of physical effects in this background could serve as a handle to deal with more complicated 
geometries like Schwarzschild. The Rindler geometry shares most of the qualitative features 
of black holes and is simple enough to allow detailed analysis. Another motivation for the 
investigation of quantum effects in the Rindler space is related to the fact that this space is 
conformally related to de Sitter space and to Robertson- Walker space with negative spatial 
curvature. As a result the expectation values of the energy-momentum tensor for conformally 
invariant fields and for corresponding conformally transformed boundaries on the de Sitter and 
Robertson- Walker backgrounds can be generated from the corresponding Rindler counterpart 
by the standard transformation (see, for instance, 

An interesting topic in the investigations of the Casimir effect is the dependence of the 
vacuum characteristics on the type of the vacuum. Vacuum expectation values of the energy- 
momentum tensor induced by an infinite plane boundary moving with uniform proper accel- 
eration through the Fulling-Rindler vacuum was studied by Candelas and Deutsch for the 
conformally coupled AD Dirichlet and Neumann massless scalar and electromagnetic fields. In 
this paper only the region of the right Rindler wedge to the right of the barrier is considered. 
In Ref. ^5] we have investigated the Wightman function and the vacuum energy-momentum 
tensor for a massive scalar field with general curvature coupling parameter, satisfying the Robin 
boundary conditions on the infinite plane in an arbitrary number of spacetime dimensions and 
for the electromagnetic field. We have considered both regions, including the one between the 
barrier and Rindler horizon. The vacuum expectation values of the energy-momentum tensor 
for scalar fields with Dirichlet and Neumann boundary conditions and for the electromagnetic 
field in the geometry of two parallel plates moving by uniform accelerations are investigated in 
Ref. [I3j- In particular, the vacuum forces acting on the boundaries are evaluated. They are 
presented as a sum of the 'interaction' and self-action parts. The 'interaction' forces between the 
plates are always attractive for both scalar and electromagnetic cases. Due to the well-known 
surface divergences in the boundary parts, the total Casimir energy cannot be obtained by di- 
rect integration of the vacuum energy density and needs an additional renormalization. In Ref. 
|14j by using the zeta function technique, the Casimir energy is evaluated for massless scalar 
fields under Dirichlet and Neumann boundary conditions, and for the electromagnetic field with 
perfect conductor boundary conditions on one and two parallel plates. On background of man- 
ifolds with boundaries, the physical quantities, in general, will receive both volume and surface 
contributions and the surface terms play an important role in various branches of physics. An 
expression for the surface energy-momentum tensor for a scalar field with general curvature 
coupling parameter in the general case of bulk and boundary geometries is derived in Ref. |15j . 
In Ref. the vacuum expectation value of the surface energy-momentum tensor is evaluated 
for a massles scalar field obeying Robin boundary condition on an infinite plane moving by 
uniform proper acceleration. By using the conformal relation between the Rindler and de Sitter 
spacetimes and the results from ^2], in Ref. the vacuum energy-momentum tensor for a 
scalar field is evaluated in de Sitter spacetime in presence of a curved brane on which the field 
obeys the Robin boundary condition with coordinate dependent coefficients. 

In the present paper the Wightman function and the vacuum expectation value of the field 
square the energy-momentum tensor are investigated for a massive scalar field with an arbitrary 
curvature coupling parameter obeying the Robin boundary conditions on two parallel branes 
moving by uniform proper accelerations through the Fulling-Rindler vacuum. The general case 
is considered when the constants in the boundary conditions are different for separate plates. 
Robin type conditions are an extension of Dirichlet and Neumann boundary conditions and 
appear in a variety of situations, including the considerations of vacuum effects for a confined 
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charged scalar field in external fields ^Hl> spinor and gauge field theories, quantum gravity 
and supergravity |191 !2Uj . Robin conditions can be made conformally invariant, while purely- 
Neumann conditions cannot. Thus, Robin-type conditions are needed when one deals with 
conformally invariant theories in the presence of boundaries and wishes to preserve this imparl- 
ance. It is interesting to note that the quantum scalar field satisfying the Robin condition on 
the boundary of cavity violates the Bekenstein's entropy-to-energy bound near certain points in 
the space of the parameter defining the boundary condition |2j. The Robin boundary condi- 
tions are an extension of those imposed on perfectly conducting boundaries and may, in some 
geometries, be useful for depicting the finite penetration of the field into the boundary with the 
'skin-depth' parameter related to the Robin coefficient. Robin boundary conditions naturally 
arise for scalar and fermion bulk fields in the Randall-Sundrum model |22[ I23 | l2^j . In this model 
the bulk geometry is a slice of anti-de Sitter space and the corresponding Robin coefficients are 
related to the curvature scale of the space. 

The outline of this paper is the following. In the next section the Wightman function is 
considered. The corresponding mode-sum is evaluated by using the generalized Abel-Plana 
summation formula |23]. This allows us to extract from the corresponding vacuum expectation 
values the Wightman function for the geometry of a single plate and to present the remained 
part in the form of the exponentially convergent integrals. The vacuum expectation values of the 
field square and the Casimir energy- momentum tensor is evaluated in Section |21 Various limiting 
cases are considered. In Section |3] we investigate the vacuum 'interaction' forces between the 
plates as functions on corresponding proper accelerations. Section [5] contains a summary of the 
work and some suggestions for further research. In Appendix El on the base of the generalized 
Abel-Plana formula, a summation formula is derived for the series over zeros of a combination 
of the Bessel modified functions with an imaginary order. 



2 Wightman function 

We consider a real scalar field <p(x) with general curvature coupling parameter £ satisfying the 
field equation 

(V^V + m 2 + (R) if = 0, (2.1) 

where R is the scalar curvature for a {D + l)-dimensional background spacetime, and is the 
covariant derivative operator. For special cases of minimally and conformally coupled scalars 
one has £ = and Q = {D — 1)/4D, respectively. Our main interest in this paper will be 
the Wightman function, the vacuum expectation values (VEVs) of the field square and the 
energy-momentum tensor in the Rindler spacetime induced by two parallel plates moving with 
uniform proper acceleration when the quantum field is prepared in the Fulling-Rindler vacuum. 
For this problem the background spacetime is flat and in Eq. (|2.1|) we have R = 0. As a 
result the eigenmodes are independent on the curvature coupling parameter. However, the 
local characteristics of the vacuum such as energy density and vacuum stresses depend on this 
parameter. 

In the accelerated frame it is convenient to introduce Rindler coordinates (r, £, x) related to 
the Minkowski ones, (t, x ,x) by formulas 

t = £sinliT, x 1 = £coshr, (2-2) 

where x = (x 2 , . . . , x D ) denotes the set of coordinates parallel to the plates. In these coordinates 
the Minkowski line element takes the form 

ds 2 = fdr 2 - df - dx 2 , (2.3) 
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and a world line denned by £, x = const describes an observer with constant proper acceleration 
The Rindler time coordinate r is proportional to the proper time along a family of uniformly 
accelerated trajectories which fill the Rindler wedge, with the proportionality constant equal to 
the acceleration. Assuming that the plates are situated in the right Rindler wedge x 1 > \t\, we 
will let the surfaces £ = a and £ = b, b > a represent the trajectories of these boundaries, which 
therefore have proper accelerations a -1 and 6 _1 (see Fig. We will consider the case of a 
scalar field satisfying Robin boundary conditions on the surfaces of the plates: 

= 0, j = a,b, (2.4) 

with constant coefficients Aj and Bj. Dirichlet and Neumann boundary conditions are obtained 
from here as special cases. All results below will depend, of course, on the ratios Aj/Bj only. 
However, to keep the transition to Dirichlet and Neumann cases transparent, we write the 
boundary conditions in the form l|2.4|) . 





t 




= const / 


\. \t = a V b 



Figure 1: The (x , t) plane with the Rindler coordinates. The heavy lines £ = a and £ = b 
represent the trajectories of the plates. 

The plates divide the right Rindler wedge into three regions: < £ < a, £ > b, and 
a < £ < b. The VEVs in two first regions are the same as those induced by single plates 
located at £ = a and £ = b, respectively. As these VEVs are investigated in Ref. ^21) i n the 
consideration below we restrict ourselves to the region between the plates. First we consider 
the positive frequency Wightman function G + (x,x') = (0 \(p(x)<p(x')\ 0), with |0) being the 
amplitude for the corresponding vacuum state. The VEVs of the field square and the energy- 
momentum tensor can be evaluated on the base of this function. In addition, the Wightman 
function determines the response of the particle detector of Unruh-DeWitt type, moving through 
the vacuum under consideration. By expanding the field operator over the complete set of 
eigenfunctions {(p a (x), <p* a {xj} satisfying boundary conditions (j2.4j) and using the commutation 
relations one finds 

G+(x,x') = J2<P*(x)v*a(x'), (2-5) 

a 

where the collective index a is a set of quantum numbers specifying the solution. 

To evaluate the mode sum in formula Q2.5JI we need the form of the eigenfunctions (f a (x) 
(for a recent discussion of eigenmodes in four Rindler sectors and relations between them see, 
for instance, 26 ). For the geometry under consideration the metric and boundary conditions 
are static and translational invariant in the hyperplane parallel to the plates. It follows from 
here that the corresponding part of the eigenfunctions can be taken in the standard plane wave 
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form: 

ip a = C<p{£) exp [i (kx — wr)] , a = (k, of), (2-6) 

with the wave vector k = (&2, • • • , fee)- The frequency cj in Eq. ()2.6j) corresponds to the 
dimensionless coordinate r and hence is dimensionless. The proper time r g and the frequency 
Ug measured by a uniformly accelerated observer with the proper acceleration g and world line 
(x 1 ) 2 — t 2 = g~ 2 are related to r and uj by the formulas r g = rjg, oj g = ojg (the features of 
the measurements for time, frequency, and length relative to a Rindler frame as compared to 
a Minkowski frame are discussed in Ref. |2Z])- The equation for 0(£) is obtained from field 
equation (|2.1|) on background of metric ()2.3)) and has the form 

eV'co + &'(0 + - A 2 e 2 ) 0(0 = o, (2.7) 

where the prime denotes a differentiation with respect to the argument of the function, 



A = ^k 2 + m 2 , (2.8) 

and k = |k|. In the region between the plates the linearly independent solutions to equation 
1)2.7)) are the Bessel modified functions Iiu(k£) and Kiu(k£). The solution satisfying boundary 
condition (|2.4() on the plate £ = b has the form 

4? ( A£, Aft) = Jff (A6)^ (AO - (A6)4, (A£) . (2.9) 
Here and below for a given function f(z) we use the notations 

f(j)( z )=A ] f(z) + ^zf'(z), j = a,b. (2.10) 

Note that function (|2.9|) is real, (X^, Xb) = Z^j UJ (X^, Xb). From the boundary condition on 
the plate £ = a we find that the possible values for u are roots to the equation 

Z iw (Ao,A6) = 0, (2.11) 

with the notation 

Z w (u,«) = fW^M - X"«(«)4 a) («). (2-12) 

For a fixed A, the equation (|2.11|) has an infinite set of real solutions with respect to uj. We will 
denote them by uj n = oj n (Xa, Xb), oj n > 0, n = 1, 2, . . ., and will assume that they are arranged 
in the ascending order oj n < u) n +\. In addition to the real zeros, in dependence of the values 
of the ratios Aj/Bj, equation ()2.11j) can have a finite set of purely imaginary solutions. The 
presence of such solutions leads to the modes with an imaginary frequency and, hence, to the 
unstable vacuum. In the consideration below we will assume the values of the coefficients in Eq. 
(j2.4|) for which the imaginary solutions are absent and the vacuum is stable. 

The coefficient C in formula (|2.6|) is determined from the standard Klein-Gordon orthonor- 
mality condition for the eigenfunctions which for metric 1)2.3 J) takes the form 

r b d£ <— > 

dx / — ip a d T (p* a , = id aa i. (2.13) 



(2.14) 



The ^-integral on the left of this formula is evaluated using the integration formula 

— <t>\u,Kt)<hv(£) = £ , , 2 _„2 



a i ' UJ 2 -V 2 
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valid for any two solutions <^(£), I = 1,2, to equation (|2.7|) . Taking into account boundary 
condition ()2.4j) . from Eq. (|2.13|) for the normalization coefficient one finds 



C 2_ i2^- D I ( £i\a) 



lS(Xb)i;Z iu (Xa,\b) 
Now substituting the eigenfunctions 

<p a (x) = CZ^jXt, Xb) exp [i (kx - u n r)} (2.16) 
into the mode sum formula (j2.5ft . for the positive frequency Wightman function one finds 



(2.15) 



_ f dkeW*-*') f> 4 a) (Aa)e—(^') 



(2^) D - 1 ^/£>(A6)^(Aa,A6) 
x^(A£,A6)z£>(A£',A6) . (2.17) 

Ul=OJ n 

As the expressions for the eigenfrequencies u n (as functions on Ao, Xb, Aj/Bj) are not explicitly 
known, the form ()2.17j) of the Wightman function is inconvenient. For the further evoluation of 
this VEV we can apply to the sum over n the summation formula (|A.4|) derived in Appendix lAl 
on the base of the generalized Abel-Plana formula. As a function F{z) in formula (|A.4|) let us 
choose 

_ z<»(\ ( ,\b)zi';>(\e,\b) lT _ Tl) 

Condition ()A.5|) for this function is satisfied if a 2 e' T ~ T ' < In particular, this is the case 
in the coincidence limit r = r' in the region under consideration: £, £' > a. By using formula 
(|A.4|) . for the Wightman function one obtains the expression 



dke ik(x-x') 

7r(2vr) J 

xZ^(A£, A6)zW(Ae', Aft) cosher - r% (2.19) 



G+^a:') = G + (x,x';b)- \ _ D _ x \ du Q^Xa, Xb) 



o 



where we have introduced the notation 



0^(Ao,A6) = " lAQJ . (2.20) 

/i b) (A6)Z w (Aa,A6) 



InEq. lt2~T31) 

dke ik(x-x') /.oo 



G + (x,x';b) = J ^2(2^)^-1 J dusinH™) 



vr -i»(T-r> ^{xn,xb)z%{xg>,xb) 

is the Wightman function in the region £ < b for a single plate at £ = 6. This function is 
investigated in Ref. ^2] an d can be presented in the form 

G + (x, x'; b) = G+(x, x') + ((p(x)<p(x')) (b) , (2.22) 
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where G~^(x,x') is the Wightman function for the right Rindler wedge without boundaries and 
the part 



vr(2^-i J j(0) (A6) 
xI w (Af)I w (AO cosher - r')] (2.23) 

is induced in the region £ < 6 by the presence of the plate at £ = 6. Note that the representation 
(|2.22[) with (|2.23f) is valid under the assumption ££' < b 2 & T ~ T I Hence, the application of the 
summation formula based on the generalized Abel-Plana formula allowed us (i) to escape the 
necessity to know the explicit expressions for the eigenfrequencies u n , (ii) to extract from the 
VEVs the purely Rindler and single plate parts, (iii) to present the remained part in terms of 
integrals with the exponential convergence in the coincidence limit. 
By using the identity 

+ nWn iw (Ao, A6)i#)(A£, \j)Z®(\?, Xj), (2.24) 

j=a,b 

with = 1, n" = — 1, the Wightman function can be also presented in the form 



G + (x,x') = G + (x,x';a)- j _ f ^_ / dw Q^Xa, Xb) 



dk e ik ( x_x ') f°° 
7r(2vr) J 

xZ^ a) (A£, Aa)zi a) (Ae', Aa) cosher - r')]. (2.25) 

In this formula 

G + (x, x'; a) = G+{x, x') + (ip(x)<p(x')) {a) (2.26) 
is the Wightman function in the region £ > a for a single plate at £ = a, and 



vr(2vr)^-i 7 /^ a) (Aa) 

x^(A£)^(A£') cosh[w(r - t% (2.27) 

In Eq. Q2.24|) we use the notations 

4 a) (A£,Aa) = J( a )(Aa)^(AO-^ a) (Aa)/ w (Ae), (2.28) 

flUXa,\b) = _„ Ku) {Xb) . (2.29) 

id, a) (Aa)Z w (Aa,A£>) 

Two representations of the Wightman function, Eqs. (|2.19|) and 1)2.25(1 . are obtained from each 
other by the replacements 

a «± b, I^^K^. (2.30) 

In the coincidence limit the second term on the right of formula (|2.19|) is finite on the plate 
£ = b and diverges on the plate at £ = a, whereas the second term on the right of Eq. (|2.25j) 
is finite on the plate £ = a and is divergent for £ = 6. Consequently, the form ()2.19j) [ Q2.25JI ] is 
convenient for the investigations of the VEVs near the plate £ = = a). Note that in the 
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formulas given above the integration over the angular part of the vector k can be done with the 
help of the formula 



/ 



(2tt)— Jo >' 



for a given function F(k), and Ji/(<z) is the Bessel function. In this section we have considered 
the positive frequency Whightman function. By the same method any other two-point function 
(Hadamard function, Feynman's Green function, etc.) can be evaluated. 



3 Casimir densities 

3.1 VEV for the field square 

In this section we will consider the VEVs of the field square and the energy-momentum tensor in 
the region between the plates. As the corresponding quantities for a single plate are investigated 
in Ref. ^2j, here we will be concentrated on the parts induced by the presence of the second 
plate. In the coincidence limit from the formulas for the Wightman function one obtains two 
equivalent forms for the VEV of the field square: 

<o|^ 2 |o> = (o R y\o R ) + (^ 

f'OC f'OO 

-A D / dkk°- 2 / duj£l jul {Xa, Xb)Z^ 2 (X£, Xj), (3.1) 



corresponding to j = a and j = b, and \0r) is the amplitude for the Fulling- Rindler vacuum 
without boundaries, 

Ad = • (3.2) 

2^-2 7 r^ i r(^i) 

In Eq. (|3.1|) the part (ip 2 )^ is induced by a single plate at £ = j when the second plate is 
absent. From (|2~23|) . (|2?27|) for this part one has [12] 

<^ 2 > (a) = ~Ad [°°dkk D - 2 [ X duIz^KZte), (3.3a) 
Jo Jo Kr{\a) 

= -A D T dkk D - 2 T du^^lliXi). (3.3b) 
Jo Jo £>(Xb) 

The last term on the right of formula ()3.1|) is finite on the plate at £ = j and diverges for the 
points on the other plate. 

Extracting the contribution from the second plate, we can write the expression ()3.1|) for the 
vacuum expectation value in the symmetric form 

<0 |^ 2 | 0> = (0 R \^ 2 \ R ) + J2 <^> (J) + {^) (ab) , (3-4) 

j=a,b 

with the 'interference' part 

„2\(a&) _ A _ ;„D-2 



(3.5) 



yy ab > = -A D / dkk D - 2 / du;I^(Xa) 
Jo Jo 

zL b)2 (X£,Xb) K 2 (XC) 



lL b \Xb)ZUXa,Xb) K^{Xa) 
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An equivalent form for this part is obtained with the replacements (|2.3U|) in the subintegrand. 
The 'interference' term (j3.5|) is finite for all values of £ in the range a < £ < b, including the 
points on the boundaries. The well-known surface divergences are contained in the single plate 
parts only. To find the corresponding asymptotic behaviour we note that for the points near 
the boundaries the main contribution into the w-integral comes from large values of u and we 
can use the uniform asymptotic expansions for the modified Bessel functions for large values of 
the order (see, for instance, |2H])- Introducing a new integration variable x = k/u and replacing 
the modified Bessel functions by their uniform asymptotic expansions, in the limit £ — > j to the 
leading order one obtains 

^ * u ?™\T\d ? (3 - 6) 



r 2\U) 

where 



k j = l- 25 Bj0 . (3.7) 

This term has different signs for Dirichlet and non-Dirichlet boundary conditions and is the 
same as that for a plate on the Minkowski bulk with |£ — j\ being the distance from the plate. 

In the limit a — * b with the fixed values of the coefficients in the boundary conditions and 
the mass, the 'interference' part (|3.5jl is divergent and for small values of b/a — 1 the main 
contribution comes from large values of u. Again, introducing an integration variable x = k/tu 
and replacing the modified Bessel functions by their uniform asymptotic expansions, to the 
leading order one obtains 



{LP) iWJo k akb e^ - 1 • (3 - 8) 

Large values of the proper accelerations for the plates correspond to the limit a, b —* 0. In 
this limit the plates are close to the Rindler horizon. From formulas ()3.3a[) . ()3.3b|) . (J3.5|) we see 
that for fixed values of the ratios a/b, both single plate and 'interference' parts behave as 
foi-D j n limit 6^0. In the limit a — > for fixed values £ and b, the left plate tends to 
the Rindler horizon for a fixed world line of the right plate. The main contribution into the 
w-integral in Eq. (|3,5|) comes from small values oj, lo < l/ln(2/Aa). Using the formulas for the 
Bessel modified functions for small arguments, it can be seen that the 'interference' part (|3.5f) 
vanishes as ln~ 2 (26/a). 

Now we turn to the limit of small accelerations of the plates: a, b — > oo with fixed values 
b — a, Bj/Aj, and m. In this case the main contribution comes from large values of u. Using the 
uniform asymptotic formulas for the Bessel modified functions, the following formula is obtained 
for the single plate parts: 

r(f) J m Cj{y) 



with the notation 



Aj + nWBjy 
Similarly, for the 'interference' term we find 



fa) = 4^fe ^ = 1, „W = -1. (3.10) 



D 



D 



r(f) J m Q,(!/)c5(!,)e 2 »(»-») - 1 



j=a,b 



-2j/|?-i| 

cj(y) 



(3.11) 



Formulae ()3.9|) and ()3.11|) coincide with the corresponding expressions for the geometry of two 
parallel plates on the Minkowski bulk. In this limit, £ corresponds to the Cartesian coordinate 
perpendicular to the plates which are located at £ = a and £ = b. 
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For large values of the mass, ma S> 1, we introduce in (|3.3a|) and (l3~3bl a new integration 
variable y = X/m. The main contribution into the w-integral comes from the values u ~ ^/ma. 
By using the uniform asymptotic expansions for the Bessel modified functions for large values 
of the order and further expanding over u/ma, for the single plate parts to the leading order 
one finds 

\f ) — ; — r° — r~, r^' 

2(47r)T Cj ( m )|^- j\2 
for j = a, b. By the similar way, for the 'interference' part we obtain the formula 



(3.12) 



2\(ab) 



m 



f-l e 2m(a-6) 



(47r)? ^c a (m)c b {m)(b - a)? 



(3.13) 



As we could expect, the both single plate and 'interference' parts are exponentially suppressed 
for large values of the mass. 

3.2 VEV of the energy- momentum tensor 

By using the field equation it can be seen that the expression for the energy-momentum tensor 
of the scalar field under consideration can be presented in the form 



and the corresponding trace is equal to 

T? = D(( - Cc)ViVV 2 + mV- 



(3.14) 



(3.15) 



By virtue of Eq. (|3.14|) . the VEV of the energy- momentum tensor is expressed in terms of the 
Wightman function as 



(0\T ik (x)\ 0) = lim ViV k G + (x,x') + 



C 



(0|(^ 2 (x)|0). (3.16) 



Making use the formulas for the Wightman function and the field square, one obtains two 
equivalent forms, corresponding to j = a and j = b (no summation over i): 



<0|7f|0> = (0 R \Jt\0 R ) + (lt) U) -A D 8* J ' dkk D ~ 2 

POD r 

xA 2 / dton ju} (Xa,Xb)F^ zW(X£,\j) 
Jo L 



In this formula, 



(0 R \T?\0 R ) = 6t-?- \ 
K Jo 



dkk°- 2 X 2 



dusmlnrujf {i) [K tuJ {X^} 



(3.17) 



(3.18) 



o 



is the corresponding VEV for the Fulling-Rindler vacuum without boundaries, and the terms 
(no summation over i) 



(T t k ) (a) 
{T k ){b ) 



-Ad 8. 



-A D 5 



dkk°- 2 X 2 



dkk°- 2 X 2 



o K™(\a) 
' K^(Xb) 

du^rT\ 

o li b) (Xb) 



F®[I u (>t)], 



(3.19a) 
(3.19b) 
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are induced by the presence of a single plane boundaries located at £ = a and £ = b in the 
regions £ > a and £ < b respectively. In formulas ()3.17j) . (|3.19aj) . ()3.19b|) for a given function 
g(z) we use the notations 



F(%(z)] 
F(%(z)] 
F®\g(z)] 



1 



2C 



+ 1 + 



or 



5 2 W 



2 cZ^ " 



+ ^ 5 2 (z)-^ 2 (z), (3.20) 



1 / rfg(z) 

2 I 



3"* 



C d 2 



dgjz) 
dz 



+ {!+ \g\ Z ) 



9 2 {z) k 2 
D - 1 A 2 ' 



(3.21) 
(3.22) 



where i = 2, . . . , D and the indices 0,1 correspond to the coordinates r, £, respectively. For the 
last term on the right of Eq. (|3.17|) we have to substitute g{z) = zjj\z, Xj). The expressions for 
the functions f^[g(z)] in Q3.18[) are obtained from the corresponding expressions for F®\g(z)] 
by the replacement to — > iu. It can be easily seen that for a conformally coupled massless scalar 
the energy-momentum tensor is traceless. 

The purely Fulling-Rindler part (|3.18l) of the energy-momentum tensor is investigated in a 
large number of papers (see, for instance, references given in ^Sl)- The most general case of 
a massive scalar field in an arbitrary number of spacetime dimensions has been considered in 
Ref. |29| for conformally and minimally coupled cases and in Ref. |12j for general values of 
the curvature coupling parameter. For a massless scalar the VEV for the Rindler part without 
boundaries can be presented in the form 



(Ti 



*\(«) 

i / sub 



(0r\T?\0r) 



(Q M \7i\0 M ) 



(4*)*r(§) 



uj D g^(uj)dw 



+ (-1 



(3.23) 



where the expressions for the functions g®(u) are presented in Ref. ^21 > an d |0m) is the 
amplitude for the Minkowski vacuum without boundaries. Expression (|3.23|) corresponds to the 
absence from the vacuum of thermal distribution with standard temperature T = (2-7r£) _1 . In 
general, the corresponding spectrum has non-Planckian form: the density of states factor is 
not proportional to LO D ~ 1 dto. The spectrum takes the Planckian form for conformally coupled 

lu) = 1, i = 1,2, ... D. It is of interest to note that 



scalars in D = 1,2,3 with g^°\to) = —Dg 



for even values of spatial dimension the distribution is Fermi-Dirac type (see also |3()| 131)). For 
the massive scalar the energy spectrum is not strictly thermal and the corresponding quantities 
do not coincide with ones for the Minkowski thermal bath. 

The boundary induced quantities (j3.19a|) . ()3.19b|) are investigated in Ref. for a confor- 
mally coupled D = 3 massless Dirichlet scalar in the region on the right from a single plate and 
in Ref. |12j for a massive scalar with general curvature coupling and Robin boundary condition 
in an arbitrary number of dimensions in both regions. The single boundary parts diverge at 
the plates surfaces £ = j, j = a,b. Near the plates the leading terms of the corresponding 
asymptotic expansions have the form (no summation over i) 



(D-l)(j-0 



2^r^|£ 



\D+1 



J ' 



(3.24) 



with i = 0, 2, . . . , D, and kj is defined by Eq. Q3.7j) . These leading terms vanish for a conformally 
coupled scalar and coincide with the corresponding quantities for a plane boundary in the 
Minkowski vacuum. 
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Now let us present the VEV (|3.17|) in the form 



(0\T z k \0) = (0 R \T t k \0 R ) + ]T (T?)^ + {7*)W, 



(3.25) 



where (no summation over i) 

^rpk^(ab) 



-A D 5 



]=a,b 



dk k°- 2 X 2 
o ./o 

(6)/ 



dwI^XXa) 



F®[Z$\K,\b)\ FW[^(AO] 



li fe) (A6)Z w (Aa,A6) K^(Aa) 



(a), 



(3.26) 



is the 'interference' term. The surface divergences are contained in the single boundary parts 
and this term is finite for all values a < £ < b. An equivalent formula for (7^}( a6 ) is obtained 
from Eq. (|3.26|) by replacements ()2.3(J|) . 

Both single plate and 'interference' parts separately satisfy the standard continuity equation 
for the energy- momentum tensor, which for the geometry under consideration takes the form 



/rpO\ 
\ J / 



(3.27) 



For a conformally coupled massless scalar field the both parts are traceless and we have an 
additional relation (T-) = 0. 

In the limit a — * b expression (|3.26|) is divergent and for small values of b/a— 1 the main 
contribution comes from the large values of u>. Introducing a new integration variable x = k/uj 
and replacing Bessel modified functions by their uniform asymptotic expansions for large values 
of the order, at the leading order one receives 



i\(ab) 



(47T) 



D 
2 



r(# 



dy 



V 



D 



k a k b e 2 y( b - a *> - 1 



1 + 2D (1 - 6\) (£ - £ c 



j=a,b 



(3.28) 



In the limit of large proper accelerations for the plates, a,b — ► 0, for fixed values a/b and 
£/&, the world lines of both plates are close to the Rindler horizon. In this case the single plate 
and 'interference' parts grow as b~ D ~ l . The situation is essentially different when the world line 
of the left plane tends to the Rindler horizon, a — > 0, whereas b and £ are fixed. By the way 
similar to that for the case of the field square, it can be seen that in this limit the 'interference' 
part (|3.26f) vanishes as ln _2 (26/a). 

In the limit of small proper accelerations, a, b — > oo with fixed values b — a, Bj/Aj, and 
m, the main contribution comes from large values of u. Using the asymptotic formulas for the 
Bessel modified functions, to the leading order one obtains (no summation over i) 



4(1-4) 
(4vr)fr(#; 

e -2y|5-il 



dy (y 2 

c- 



2\ — 

m 2 



C-Cc + 



1/4 



(3.29) 
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for the single boundary terms, and 



Pi 



i\ (ab) 



(47T) 



— s- roc 



2 + 1 ) 
" (1 " *l) 



dy- 



(y 2 — m 2 ) 1 



c a {y)c b {y)e 2 y( b - a ) 

4£>(£ - C c )y 2 - m 2 



2(y 2 — m 2 ) 



cj(y) 



(3.30) 



for the 'interference' term and with the function Cj(y) defined by (j3.10j) . These expressions are 
exactly the same as the corresponding expressions for the geometry of two parallel plates on the 
Minkowski background investigated in [33] for a massless scalar and in Ref. [33] for the massive 
case. In particular, the single boundary terms vanish for a conformally coupled massless scalar. 

In the large mass limit, ma S> 1, by the method similar to that used in the previous sub- 
section for the field square, it can be seen that the both single plate and 'interference' parts are 
exponentially suppressed (no summation over i): (T-)^ ~ m D / 2+1 exp[— 2m|£ — j = a,b, for 



single plate parts and (T { 



i\(ab) 



~ m 



D / 2+l exp[2m(a — b)] for the 'interference' part. 



4 'Interaction' forces between the plates 

Now we turn to the 'interaction' forces between the plates due to the vacuum fluctuations. The 
vacuum force acting per unit surface of the plate at £ = j is determined by the {-component 
of the vacuum energy-momentum tensor evaluated at this point. The corresponding effective 
pressures can be presented as a sum of two terms: 

P U) =P? +4l t) , j = a,b. (4.1) 

The first term on the right is the pressure for a single plate at £ = j when the second plate 
is absent. This term is divergent due to the surface divergences in the subtracted vacuum 
expectation values and needs additional renormalization. This can be done, for example, by 
applying the generalized zeta function technique to the corresponding mode sum. This procedure 
is similar to that used in Ref. ^3] for the evaluation of the total Casimir energy in the cases of 
Dirichlet and Neumann boundary conditions and in Ref. j!6j for the evaluation of the surface 
energy for a single Robin plate. This calculation lies on the same line with the evaluation of the 
total Casimir energy and surface densities and will be presented in the forthcoming paper |32j . 
Note that in the formulae for the VEV of the energy-momentum tensor the Robin coefficients 
enter in the form of the dimensionless combination (5j = Bj/(jAj). As a result in the massless 
case from the dimensional arguments we expect that the single plate part will have the form 
Pi = a((3j)j~( D+l \ The coefficient a(f3j) in this formula will change if we will change the 
renormalization scale and can be fixed by imposing suitable renormalization conditions which 
relates it to observables. 

The second term on the right of Eq. (|4.1j) , 

(4.2) 

with j, I = a, b, I j, is the pressure induced by the presence of the second plate, and can be 
termed as an 'interaction' force. This term is finite for all nonzero distances between the plates 
and is not affected by the renormalization procedure. Note that the term 'interaction' here 
should be understood conditionally. The quantity determines the force by which the scalar 
vacuum acts on the plate due to the modification of the spectrum for the zero-point fluctuations 
by the presence of the second plate. As the vacuum properties are <!; -dependent, there is no a 



0) 

P(int) 
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priori reason for the 'interaction' terms (and also for the total pressures pW) to be equal for 
j = a and j = b, and the corresponding forces in general are different. For the plate at £ = j 
the 'interaction' term is due to the third summand on the right of Eq. (|3.17|) . Substituting into 
this term £ = j and using the Wronskian for the modified Bessel functions one has 



Jj) 

^(int) 



A 



D A] 



2f 



dk k 



D-2 



dto [ (X 2 j 2 + to 2 ) 0] + 4(pj - 1] Q ju> (Xa,Xb), (4.3) 



with f3j defined in the paragraph after formula (|4.1j) . In dependence of the values for the 
coefficients in the boundary conditions, the effective pressures (|4,3|) can be either positive or 
negative, leading to repulsive or attractive forces. It can be seen that for Dirichlet boundary 



,Cj) 



> and 



condition on one plate and Neumann boundary condition on the other one has pY{ nt ) 
the 'interaction' forces are repulsive for all distances between the plates. Note that for Dirichlet 
or Neumann boundary conditions on both plates the 'interaction' forces are always attractive 
[T3] . By using the relation 







[(X 2 j 2 + u 2 ) (3 2 + - 1] ApjUXa, Xb) = n^j— In 



ti a) (Xa)KL b) (Xb) 



I {b) {Xb)K^{Xa) 



(a). 



(4.4) 



with rfi) from l|3~Tfi|) . expressions (|4,3|) for the 'interaction' forces can be written in another 
equivalent form 



Jj) 
P(int) 



n 



2j 



9 , 



dkk°- 2 / dto 
Jo 

lL a) (Xa)KL b) (Xb) 



1 + 



(4( - 1) Pj 



(X^f + ^)f3 2 + P 3 -l 



lL b) (Xb)K^(Xa) 



(a). 



(4.5) 



For Dirichlet and Neumann scalars the second term in the square brackets is zero. To clarify the 
dependence of the vacuum 'interaction' forces on the parameters a, b it is useful to write down 
the corresponding derivatives: 



n 



dl 



ApA 2 a A 2 
Zabj jo 
(A 2 Z 2 + a; 2 )A 2 + ft 
Z 2 (Xa,Xb) 



poo poo 

/ dk k D - 2 / dco [(X 2 j 2 + u; 2 ) (3 2 + 469,- - l] 
Jo Jo 



(4.6) 



with j,l = a,b, j / I. 

Now we consider the limiting cases for the 'interaction' forces between the plates. For small 

distances between the plates, b/a — 1 <C 1, to the leading order over 1/(6 — a), the 'interaction' 

forces are the same as for the plates in the Minkowski bulk with the distance b — a. The latter 

are determined by i = 1 component of tensor (J3.3U|) . In this limit the 'interaction' forces are 

repulsive in the case of Dirichlet boundary condition on one plate and non-Dirichlet boundary 

condition on the another, and are attractive for all other cases. Note that in the limit b — > a 

with fixed values of the boundary coefficients and the proper acceleration of the left plate, a -1 , 

(i) 

the renormalized single plate parts pY remain finite while the 'interaction' part goes to infinity. 
This means that for sufficiently small distances between the plates the 'interaction' term on the 
right of formula (|4.1j) will dominate. 

For large distances between the plates one has a/b <C 1. Introducing a new integration 
variable y = Xb and using the asymptotic formulas for the Bessel modified functions for small 
values of the argument, we can see that the subintegrand is proportional to (ya/b) 2w . It follows 
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from here that the main contribution into the u;-integral comes from small values of oj. Expanding 
with respect to w, in the leading order we obtain 



(a) 
V) 



P(int) 



ir 2 A D (1 - ACPa) A 2 
2A(D -l)a 2 b D - 1 ln 3 (2b/a 

n 2 A D A 2 



dy (y 



2,2^ y A 

m b ) 2 



2 f 



1 



nib 



2Ab D + l \n 2 (2b/c 



dyyiy 



M 
i 



T)l6 



4 6)2 ( y ) 



(4.7a) 
(4.7b) 



For a massless minimally coupled scalar field these pressures have the same sign. In Figure [21 we 
have plotted the dependence of the vacuum 'interaction' forces between the plates as functions 
on the ratio a/b for a massless minimally coupled scalar field in D = 3 with Robin coefficients 
(3 a = and /?& = 1/5. These forces are repulsive for small distances and are attractive for 
large distances. In the presented example there are parameter choices which give vanishing 
'interaction' forces. 



0.02 i 

0.015 " y<m 




0.3 



0.4 



0.5 



0.6 



Figure 2: The vacuum effective pressures aD+1 P^l t y 3 = a > b, determining the 'interaction' forces 
between the plates, as functions on a/b for a massless minimally coupled scalar field in D = 3 
with Robin coefficients (3 a = and (3b = 1/5. 

For large values of the mass, ma 3> 1, the main contribution into the w-integral comes from 
the values lo ~ ^m(b — a). By using the uniform asymptotic expansions for the Bessel modified 
functions, to the leading order one finds 

m B 2 m f+3 e 2mM)^ 

(mt) (47r)Tq(m)[A i - nO)B jm ] 2 j(b - o)f ' 

for Bj ^ 0, j, I = a,b, j ^ I. For Bj = the leading term is the same as that for the Aj = 
case. The latter is obtained from (|4.8|) . 

5 Conclusion 

The use of general coordinate transformations in quantum field theory in flat spacetime leads 
to an infinite number of unitary inequivalent representations of the commutation relations with 
different vacuum states. In particular, the vacuum state for a uniformly accelerated observer, 
the Fulling-Rindler vacuum, turns out to be inequivalent to that for an inertial observer, the 
Minkowski vacuum. In the present paper we have considered the positive frequency Wightman 



15 



function, the VEVs of the field square and the energy-momentum tensor for a scalar field in the 
region between two infinite parallel plates moving by uniform proper accelerations, assuming 
that the field is prepared in the Fulling-Rindler vacuum state and satisfies the Robin boundary 
conditions on the plates. The general case is investigated when the constants in the Robin 
boundary conditions are different for separate plates. The boundaries and boundary conditions 
are static in the Rindler coordinates and no Rindler quanta are created. The only effect of 
the imposition of boundary conditions on a quantum field is the vacuum polarization. The 
Wightman function is presented in the form of the mode sum involving series over zeros to = uj n 
of the function Zj W {Xa, Xb) defined by relation (|2.12|) . For the summation of these series we 
have applied a summation formula derived in Appendix lAl by using the generalized Abel-Plana 
formula. This allowed to extract from the Whightman function the part due to a single plate and 
to present the additional part in terms of integrals, exponentially convergent in the coincidence 
limit. The single plate part is investigated previously in Ref ^2]. The contribution induced by 
the second boundary is presented in two alternative forms, Eqs. (|2.19j) . (|2.25j) . obtained from 
each other by replacements (j2.3U|) . In Section |31 by using the expression for the Wightman 
function, we evaluate the VEVs of the field square and the energy-momentum tensor. The 
latter is diagonal and the corresponding components are determined by relation (|3.17|) . Various 
limiting cases are studied. In the limit of small distances between the plates, to the leading 
order, the VEVs are the same as those for two parallel plates in the Minkowski vacuum. In the 
near horizon limit, a,b — > 0, the proper accelerations of the plates are large. For fixed values 
a/b and £/&, the VEVs grow as b^ D for the field square and as for the components of 

the energy-momentum tensor. In the limit when the world line of the left plate tends to the 
Rindler horizon, a —* 0, for a fixed proper accelerations of the right plate and the observer, 
the VEVs induced by the left plate vanish as ln _2 (26/a) for both the field square and energy- 
momentum tensor. For large values of the mass, the both single plate and interference parts of 
the VEVs are exponentially suppressed. The vacuum forces acting on boundaries are determined 
by ^-component of the stress and are investigated in Section 0] These forces are presented as 
the sums of two terms. The first ones correspond to the forces acting on a single boundary 
then the second boundary is absent. Due to the surface divergences in the VEVs of the energy- 
momentum tensor, these forces are infinite and need an additional renormalization. The another 
terms in the vacuum forces are finite and are induced by the presence of the second boundary. 
They correspond to the 'interaction' forces between the plates. These forces per unit surface 
are determined by formula (|4.2|) . For small distances between the plates, to the leading order 
the standard Casimir result on background of the Minkowski vacuum is rederived. In this limit 
the 'interaction' forces are repulsive in the case of Dirichlet boundary condition on one plate 
and non-Dirichlet boundary condition on the another, and are attractive for all other cases. For 
large distances, the 'interaction' forces can be either attractive or repulsive in dependence of the 
coefficients in the boundary conditions. In Figure we have presented an example when the 
vacuum 'interaction' forces are repulsive for small distances and are attractive for large distances. 
This provides a possibility for the stabilization of the interplate distance by vacuum forces. 

However, it should be noted that to make reliable predictions regarding quantum stabilization, 

(i) 

the renormalized single plate parts p± also should be taken into account. The calculation of 
these quantities lies on the same line with the evaluation of the total Casimir energy and surface 
densities and will be presented in the forthcoming paper |32| . 

In the present paper we have investigated the VEV of the bulk energy-momentum tensor. For 
scalar fields with general curvature coupling and Robin boundary conditions, in Ref. j3B] it has 
been shown that in the discussion of the relation between the mode sum energy, evaluated as the 
sum of the zero-point energies for each normal mode of frequency, and the volume integral of the 
renormalized energy density for the Robin parallel plates geometry it is necessary to include in 
the energy a surface term concentrated on the boundary (see also the discussion in Ref. [35[ I36| ) . 
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Similar issues for the spherical and cylindrical boundary geometries and in braneworld scenarios 
are discussed in Refs. |371 1381 OH] . An expression for the surface energy- momentum tensor for a 
scalar field with a general curvature coupling parameter in the general case of bulk and boundary 
geometries is derived in Ref. |15j . The investigation of the total Casimir energy, the surface 
densities, and the energy balance for the geometry under consideration will be reported in |32| . 

The formulas derived in this paper can be used to generate the vacuum densities for a 
conformally coupled massless scalar field in de Sitter spacetime in presence of two curved branes 
on which the field obeys the Robin boundary conditions with coordinate dependent coefficients. 
The corresponding procedure is similar to that realized in Ref. ^Zj for the geometry of a single 
brane and is based on the conformal relation between the Rindler and de Sitter line elements. 
The results obtained above can be also applied to the geometry of two parallel plates near the 
D = 3 "Rindler wall." This wall is described by the static plane-symmetric distribution of the 
matter with the diagonal energy-momentum tensor = diag(e m , — p m , — p m , — p m ) (see Ref. 
|4()j). Below we will denote by x the coordinate perpendicular to the wall and will assume that 
the plane x = is at the center of the wall. If the plane x = x s is the boundary of the wall, 
when the external [x > x s ) line element with the time coordinate t can be transformed into 
form ((23)) with 

C(x) = x-x s + - , r = 27Ta s ^g 00 (x s )t. (5.1) 

ZTTt7 s 

In this formula the parameter a s is the mass per unit surface of the wall and is determined by 
the distribution of the matter: 



a s = 2 (e m + 3p ri 
Jo 



i \ i i/ 2 

*. (5.2) 



For the "Rindler wall" one has g' 22 (x)\x=o < [Ml (the external solution for the case 322( x )U=o > 
is described by the standard Taub metric). Hence, the Wightman function, the VEVs for 
the field square and the energy-momentum tensor in the region between two plates located at 
x = X\ and x = x%, x$ > x s near the "Rindler wall" are obtained from the results given above 
substituting £j = i = 1, 2 and £ = £(x). For a s > 0, x > x s one has > £(x s ) > and 

the Rindler metric is regular everywhere in the external region. 
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A Summation formula over zeros of Z iz (u,v) 

As we have seen in section [21 the Whightman function for a scalar field in the region between 
two plates uniformly accelerated through the Fulling-Rindler vacuum is expressed in terms of 
series over the zeros z = lo u of the function Zi z (u, v ), v > u, defined by formula (|2.12|) . To derive 
a summation formula for these series, we choose in the generalized Abel-Plana formula j^Sj the 
functions 

2i 

f(z) = — sinh-7rz F(z), 

9{z) = — r F(z), 

Zi z {lL,V) 
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with a meromorphic function F{z) having poles z = 2% (7^ uj n ), ImZf. 7^ 0, in the right half-plane 
Re z > 0. The zeros oj n are simple poles of the function g(z). By taking into account the relation 

f(^) ± /(*) = % , v F(z), (A.2) 

Z iz (u,v) 



for the function R[f(z),g(z)] in the generalized Abel-Plana formula one obtains 
R[f(z),g(z)} = 2vri 



' 00 7 (fe) (v)I {a) (u) 

?1=1 Wz^iziu, v)\z=U) n 



(A.3) 



where the zeros w n are arranged in ascending order, and = sgn(Imz ( ! c ). Substituting 

relations (|A.2|) and (jA.3|) into the generalized Abel-Plana formula, as a special case the following 
summation formula is obtained: 



d 



1 m z **( u > v )\* 



1 f°° 

F{oj n ) = — - I dz sinh irz F(z) 
Jo 



VRcg F{Z) I {b) (v)I (a) 



k 

00 



z=z k Z iz (u,v) Mz k )z^ > -ia(z k )z 







/''"'Y'f' fflMgM. (A.4) 



Here the condition for the function F(z) is easily obtained from the corresponding condition in 
the generalized Abel-Plana formula by using the asymptotic formulas for the Bessel modified 
function and has the form 



\F{z)\ <e{\z\)e- nRez (-) 2llmZl , Rez>0, (A.5) 



for \z\ — ► 00, where |z|e(|z|) — ► when \z\ — > 00. Formula (|A.4|) can be generalized for the 
case when the function F{z) has real poles, under the assumption that the first integral on the 
right of this formula converges in the sense of the principal value. In this case the first integral 
on the right of Eq. (|A.4|) is understood in the sense of the principal value and residue terms 
from real poles in the form Res<?(.z), Imz^. = 0, have to be added to the right-hand side of 

z= z k 

this formula, with g(z) from (|A.1|) . For the case of series in Eq. (|2.17|) the function F(z) is an 
analytic function and the residue terms in Eq. (|A,4|) are absent. 
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